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Spectral properties of quarks above T c in quenched lattice QCD 
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We analyze the quark spectral function above the critical temperature for deconfinement in 
quenched lattice QCD using clover improved Wilson fermions in Landau gauge. We show that 
the temporal quark correlator is well reproduced by a two-pole approximation for the spectral func- 
tion and analyze the bare quark mass dependence of both poles as well as their residues. In the 
chiral limit we find that the quark spectral function has two collective modes which correspond to 
the normal and plasmino excitations. At large values of the bare quark mass the spectral function 
is dominated by a single pole. 
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To explore the properties of hot and dense matter 
formed by quarks and gluons above the critical tempera- 
ture for deconfinement (T c ) is an intriguing problem that 
has been addressed in many studies. Recent experimen- 
tal results on the properties of the matter created in 
heavy ion collisions at the Relativistic Heavy-Ion Col- 
lider (RHIC) suggest, that its time evolution above T c 
is well described by ideal hydrodynamics down to the 
freeze-out temperature in the vicinity of T c [lj . In order 
to understand better the structure of matter in this non- 
perturbative region, it is desirable to identify the basic 
degrees of freedom of the system and their quasi-particle 
properties. 

At asymptotically high temperatures almost free 
quarks and gluons are most certainly the basic degrees of 
freedom that control the properties of the Quark-Gluon 
plasma (QGP) [2J. In this regime properties of the QGP 
can be analyzed using perturbative techniques. At lower 
temperatures the application of hard-thermal loop (HTL) 
resummation Q still allowed to define gauge independent 
propagators for quarks and gluons that can be used to 
study properties of the QGP perturbatively. From these 
perturbative analyses it is known that the collective exci- 
tations of gluons and quarks develop a mass gap (thermal 
mass) that is proportional to gT 0, H 0] j where g and T 
denote the gauge coupling and temperature, respectively. 
Moreover, the number of poles in the finite temperature 
quasi-particle propagators is doubled. In addition to the 
normal modes, which reduce to poles in the free particle 
propagator, plasmon and plasmino modes appear. 

At temperatures in the vicinity of T c it is apriori not 
clear whether a quasi-particle picture for quarks and glu- 
ons is valid at all. However, lattice results on e.g. baryon 
number and electric charge fluctuations in the vicinity 
of T c [5( suggest that quasi-particles with quark degrees 
of freedom are the carriers of these quantum numbers. 
Quasi-particles also have been used successfully to de- 
scribe lattice QCD results on the equation of state 
Moreover, the apparent quark number scaling of the el- 
liptic flow observed in the RHIC experiments [7[ may also 



suggest that quasi-particles with quark quantum num- 
bers exist even close to T c . Despite the problem of gauge 
dependence of quark and gluon propagators, it therefore 
is desirable to analyze their properties at high tempera- 
ture through a direct calculation within the framework 
of QCD. In this Letter, we analyze dynamical proper- 
ties of quarks above T c in quenched lattice QCD. These 
calculations have been performed in Landau gauge. So 
far, there have been only a few studies that address this 
problem in lattice calculations [8|. 

In order to understand the origin of the plasmino mode 
in the quark propagator in the high temperature limit, it 
is instructive to consider the quark propagator at inter- 
mediate temperature by introdu cing some energy scale of 
the order of the temperature 

BUSED]- In 0> the tem- 
perature dependence of the spectral function for fermions 
with scalar mass m has been considered in QED. In this 
model, the spectral function at zero temperature has two 
poles at energies ui = ±ra, while in the high temperature 
limit, T/m — > oo, it approaches the HTL result, which 
has four poles. The one-loop calculation performed in Q 
clearly showed that the two limiting forms of the spec- 
tral function are connected continuously; in the spectral 
function a peak corresponding to the plasmino gradually 
appears and becomes larger with increasing tempera- 
ture, in addition to the normal quasi-particle peak [9(. 

In this Letter, we analyze the quark propagator at two 
values of the temperature, T = 1.5T C and 3T C , as a func- 
tion of the bare quark mass. To simplify the present 
analysis, all our calculations have been performed for zero 
momentum. The dynamical properties of quarks at zero 
momentum are encoded in the quark spectral function 
p(tjj) which is related to the Euclidean correlation func- 
tion 
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d 3 xd 3 y{iP(r,x)^(0,y)), 



through an integral equation 
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3 16 64,48 7.457 1.3389 0.13390 0.015 

12 48 7.192 1.3550 0.13437 0.021 

1.5 16 64,48 6.872 1.4125 0.13495 0.031 

12 48 6.640 1.4579 0.13536 0.041 



TABLE I: Simulation parameters [If 



with the quark field ■0, the spatial volume V, and the 
imaginary time r which is restricted to the interval < 
r < 1/T. The Dirac structure of p{ui) is decomposed as 



= p + HA + 7° + p_HA_7 , 



(3) 



with projection operators A± = (1 ± 7°)/2. The charge 
conjugation symmetry leads to Pq{uj) = po(—u>), p & 
-p„(-u), and p+(u>) = p~(-u) = p (uj) + p s {uj) 
In the following analysis, we concentrate on a determi- 
nation of p±(u>) instead of po tB (u>), because excitation 
properties of quarks are more apparent in these chan- 
nels [9(. In analogy to Eq. ([3]) we introduce the de- 
composition of the correlation function S(t) as S(t) — 
5+(r)A + 7° + S , _(t)A_7°, where S± are related through 
S+{t)=S-(P-t). 

For free quarks with scalar mass m q the spectral func- 
tions, p±(u>) = S(uj T Tn q ), have quark and anti-quark 
poles at uj = ±to 9 , respectively. In the high temperature 
limit, additional poles, corresponding to the plasmino, 
appear at negative energy for p+(u>) and positive energy 
for p-(oS) f^. Txxl. [l^| . While the positivity of p±{oS) is en- 
sured by definition, these spectral functions are neither 
even nor odd functions. In the chiral limit, however, p s 
vanishes and p±(u>) become even functions. 

To extract the spectral function p+(<jj) from S(t) using 
Eq. (|2|), we assume that p+(oj) can be described by a two- 
pole ansatz, 



p+(w) = Zi6(lu - Ei) + Z 2 5(uj + E 2 ), 



(4) 



where the residues and energies Ei, 2 > have to 
be determined from a fit to S+(t). The poles at lo = 
Ei, —Ei correspond to the normal and plasmino modes, 
respectively 

The correlation function S(t) has been calculated at 
two values of the temperature, T = 1.5T C and 3T C , in 
quenched QCD using non-perturbatively improved clover 
Wilson fermions [3, To control the dependence 

of our results on the finite lattice volume, N% x N T , 
and lattice spacing, a, we analyze the quark propaga- 
tor on lattices of three different sizes. The gauge field 
ensembles used for this analysis have been generated and 
used previously by the Bielefeld group to study screening 
masses and spectral functions [151 ]. The different simu- 
lation parameters are summarized in Table [I] |l5] . For 



each lattice size, 51 configurations have been analyzed. 
On the 64 3 x 16 lattices and at our smallest tempera- 
ture, T — 1.5T C , we observed for the largest values of 
the hopping parameter, i.e. closest to n c , an anomalous 
behavior of the quark propagator on a few gauge field 
configurations. The appearance of such exceptional con- 
figurations in calculations with light quarks in quenched 
QCD is a well-known problem in calulations with Wil- 
son fermions [lf|. We identified 7 such configurations, 
which we excluded from our analysis. The properties 
of the quark propagator on these configurations will be 
discussed in more detail elsewhere [13] • Quark propa- 
gators have been calculated after fixing each gauge field 
configuration to Landau gauge. For this we used a con- 
ventional minimization algorithm with a stopping crite- 
rion, (l/3)tr|9 M A^| 2 < 10 -11 . In the Wilson fcrmion 
formulation the bare mass, Too, is related to the hopping 
parameter k, through the standard relation 
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where k c denotes the critical hopping parameter corre- 
sponding to the chiral limit, or vanishing quark mass. 

To evaluate Eq. |T]) numerically, we solve the linear 
equation ^result = V'sourcc for a given source ^source, 
with K being the fermion matrix. For this procedure, we 
use the wall source, 0^ urce = (V^) Ei^(^ ;K )' which 
we found to be very efficient in reducing the statistical 
error in the propagator calculation. To reduce the sta- 
tistical error further, we define the correlation function 
S+(t) for each configuration by 

S+ att - (r) = ^tr[S(T)7°A + + S{£ - r) 7 °A_], (6) 

with the trace taken over Dirac and color indices. 

In Fig.[Tl we show the numerical results for S+^^t) for 
several values of k calculated on a lattice of size 64 3 x 16 at 
T = 3T C . One sees that the shape of S+ u '(t) approaches 
that of a single exponential function for smaller k, while 
it becomes symmetric as k approaches k c . In the vicinity 
of the wall source, i.e. at small and large r, we see devia- 
tions from this generic picture which can be attributed to 
distortion effects arising from the presence of the source. 
We thus exclude points with r < r m i„ and N T — r < r min 
from our fits to the ansatz given in Eq. (0| . The resulting 
correlation functions obtained from correlated fits with 
T m in = 3 are shown in Fig. [TJ One sees that S , ^- tt - is 
well reproduced by our fitting ansatz the % 2 /dof 

of our fits is between 2 and 3 at 0.1335 < n < 0.134, 
while it gradually increases as k becomes smaller than 
K = 0.1335. A similar behavior is also observed for our 
other lattice sizes 



17] 



In Fig. [21 we show the dependence of Ei y 2 and Z 2 /(Zi + 
Z 2 ) on the bare quark mass too for T = 1.5T C and 3T C . 
The results have been obtained from two-pole fits on lat- 
tices of size 64 3 x 16. Errorbars have been estimated 
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FIG. 1: The lattice correlation function S+ u -(t) at T = 3T C 
for the lattice of size 64 3 x 16 with various values of k, and 
the fitting result with the ansatz Eq. Q. 
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FIG. 3: The bare quark mass dependence of parameters E\, 



E 2 at T = 3T C 
48 3 x 12. 



for lattices of size 64 3 x 16, 48 d x 16 and 
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The bare quark mass dependence of fitting parame- 
and Z 2 /(Zi + Z 2 ) at T = 1.5T C and 3T C for lattice 



from a Jackknife analysis. The dotted line in this figure 
denotes the pole mass determined from the bare lattice 
mass given in Eq. ([5]), i.e. E\/T = ttlq/T = maaN T . 
The figure shows that the ratio Z 2 /(Zi + Z-x) becomes 
larger with decreasing too and eventually reaches 0.5. 
The hopping parameters satisfying Z\ = Z 2 are k' c — 
0.133974(10) for T = 3T C and k' c = 0.134991(9) for 
T = 1.5T C , which are consistent with the values for Kr 
given in Table HI The latter had been obtained in [15[ 
from a fit to critical hopping parameters determined in 
14j from the vanishing of the isovector axial current. The 
numerical results obtained on 64 3 x 16 lattices show that 



Ei and E 2 are equal within statistical errors at k — k' c . 
The spectral function p+(w) thus becomes an even func- 
tion at this point; the quark propagator becomes chi- 
rally symmetric despite the presence of a thermal mass, 
tot = E\ — Ei. From Fig. [2 one also finds that the 
ratio ttit/T is insensitive to T in the temperature range 
analyzed in this work, while it is slightly larger for lower 
T. 

As too becomes larger, Z 2 /(Zi + Z 2 ) decreases and 
p+{to) is eventually dominated by a single-pole. One sees 
that Ei has a minimum at mo > 0, while E 2 is an in- 
creasing function of too. In the one- loop approximation, 
the peak in p+(u>) corresponding to Ei (E 2 ) is monoton- 
ically increasing (decreasing) function of too/T 0, 17|. 
The quark mass dependence of poles found here thus is 
qualitatively different from the perturbative result. We 
find, however, that slope of E 2 as function of ttiq/T de- 
creases with increasing T. This may suggest that the 
perturbative behavior could eventually be recovered at 
much larger temperatures. 

In order to check the dependence of our results on the 
lattice spacing and finite volume, we analyzed the quark 
propagator at T — 3T C for three different lattice sizes. 
Results for E\ and E 2 are shown in Fig. [3] Compar- 
ing the results obtained on lattices with different lattice 
cut-off, a, but same physical volume, i.e. 64 3 x 16 and 
48 3 x 12, one sees that any possible cut-off dependence is 
statistically not significant in our analysis. On the other 
hand we find a clear dependence of the quark energy lev- 
els on the spatial volume; when comparing lattices with 
aspect ratio N a /N T = 3 and 4 we find that the energy 
levels, Ei >2 , drop significantly. A similar behavior is ob- 
served also at T = 1.5T C . 

The presence of a strong volume dependence of the 
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quark propagator is not unexpected. In fact, the ther- 
mal quark mass arises as collective effect of low mo- 
mentum gluons; gluons at the soft scale p < gT play 
a crucial role to give rise to the thermal mass at high 
temperatures @, Q. However, on lattices with given 
aspect ratio N a /N T low momentum gluons are cut-off. 
The lowest non- vanishing gluon momentum is, p m in/T — 
2-k(N t /N <j ) 1 which still is larger than unity on lattices 
with aspect ratio N a /N T = 4. The situation may, 
nonetheless, be somewhat better in the temperature 
range explored here as the temperature dependent cou- 
pling g(T) is larger than unity. An analysis of quark 
spectral functions on lattices with even larger spatial 
volume is needed in the future to properly control ef- 
fects of small momenta. We attempted to estimate the 
thermal mass in the V — > oo limit by extrapolating 
the results obtained for two different volumina. Defin- 
ing m T = {Z\Ei + Z 2 E 2 )/{Zi + Z 2 )| K=Kc and assum- 
ing the volume dependence of tot as rriT(N T /N a ) = 
m T (0)exp(cA^/AT3) ; we obtain m T (0)/T = 0.771(18) 
for T = 3T C and m T (0)/T = 0.800(15) for T = 1.5T C . 
This suggests that finite volume effects may still be of 
the order of 15% in our current analysis of mx/T. De- 
spite these problems, our result clearly shows that light 
quarks near but above T c have a mass gap that is of col- 
lective nature similar to that in the perturbative regime. 

In this Letter, we analyzed the quark spectral function 
at zero momentum for T = 1.5T C and 3T C as functions 
of bare quark mass too in quenched lattice QCD with 
Landau gauge fixing. We found that the two-pole ap- 
proximation for p+(uj) well reproduces the behavior of 
the lattice correlation function. It is argued that the chi- 
ral symmetry of the quark propagator is restored at the 
critical value of k and the shape of the spectral function 
at this point takes a similar form as in the high tem- 
perature limit having normal and plasmino modes with 
thermal mass my. As too is increased, p+(uj) approaches 
a single-pole structure as one can naturally deduce in- 
tuitively. The non-perturbative nature of thermal gauge 
fields is reflected in the behavior of poles as functions 
of too, which is qualitatively different from the pertur- 
bative result We also note that the ratio ttit/T de- 
creases slightly with increasing T, which is expected to 
happen at high temperature where m/T should be pro- 
portional to a running coupling g(T). Although results 
on the quark propagator are gauge dependent, we expect 
that our results for its poles suffer less from gauge depen- 
dence, because the success of the pole approximation for 
p+(u>) indicates that the quark propagator has dynami- 
cal poles near the real axis, which are gauge independent 
quantities d, [ijj]. 

In the present study, we analyzed the quark spectral 
function in the quenched approximation. Although this 
approximation includes the leading contribution in the 
high temperature limit and thus is valid at sufficiently 
high T, the validity of this approximation near T c is non- 



trivial. For example, screening of gluons due to the po- 
larization of the vacuum with virtual quark antiquark 
pairs is neglected in this approximation. The coupling 



Til 2oj ] , which may cause 



to possible mesonic excitations 
interesting effects in the spectral properties of the quark 
are not incorporated, either. The comparison of the 
quark propagator between quenched and full lattice sim- 
ulations would tell us the strength of these effects near 
T c . 

In the future it will also be interesting to use results 
on the non-perturbative structure of quark propagators 
as input for phenomenological studies of the QGP phase. 
For example, thermal properties of the charm quark 17 1 
should be useful for the understanding of properties of 
charmonia above T c [l^] . The thermal mass of light 
quarks can also be used to evaluate details of their dy- 
namics [ill ]. 

Although in this letter we limited our analysis to zero 
momentum, a proper analysis at finite momentum is 
needed for the understanding of the entire quark spectral 
function. In particular, the confirmation of the existence 
of a minimum in the plasmino dispersion relation at fi- 
nite momentum Q clearly is a challenging problem. The 
exploration of the gluon propagator is also an important 
subject of further studies. To clarify the origin of the 
quark mass dependence of E\ and E 2l which is qualita- 
tively different from the perturbative result, as well as 
the T dependence of the thermal mass, are open ques- 
tions for further numerical and analytic studies. 
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